Nonlocal correlations of the local density of states in disordered quantum Hall systems 
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Motivated by recent high-resolution scanning tunneling microscopy (STM) experiments in the 
quantum Hall regime both on massive two-dimensional electron gas and on graphene, we consider 
theoretically the disorder averaged nonlocal correlations of the local density of states (LDoS) for 
electrons moving in a smooth disordered potential in the presence of a high magnetic field. The 
intersection of two quantum cyclotron rings around the two different positions of the STM tip, 
correlated by the local disorder, provides peaks in the spatial dispersion of the LDoS-LDoS corre- 
lations when the intertip distance matches the sum of the two quantum Larmor radii. The energy 
dependence displays also complex behavior: for the local LDoS-LDoS average (i.e. at coinciding 
tip positions), sharp positive correlations are obtained for tip voltages near Landau level, and weak 
anticorrelations otherwise. 



I. INTRODUCTION 



Quantum Hall systems offer a surprising dichotomy be- 
tween very universal macroscopic properties, such as the 
near perfect quantization of the Hall conductance, and 
sample-dependent physics dominated by local imperfec- 
tions, as recently observed in several local scanning tun- 
neling spectroscopy (STS) experiments both on massive 
two-dimensional electron gasi and on grapheno^. It is, 
however, known that some degree of universality can be 
recovered by performing sample (or disorder) average of 
local quantities, and indeed theoretical predictions for 
the averaged STS local density of states (LDoS) lead to 
Gaussian behavior near the Landau levels, with an en- 
ergy width and a lineshape that depend on the width 
and correlation length of the disorder distribution, re- 
spectively^. Because the information extracted from 
transport experiments is limited, correlations of current 
(i.e., noise measurements) have been previously exam- 
ined^, successfully demonstrating the existence of frac- 
tionally charged quasiparticles for the fractional quan- 
tum Hall effect. The question we wish to raise here is 
the nature of the disorder averaged nonlocal correlations 
of local physical quantities (such as the LDoS) in the 
quantum Hall regime. Such study can, in principle, be 
experimentally achieved by sampling large spatial areas 
of the sample surface using the displacements of the STM 
tip and correlating the measured LDoS at two different 
tip positions (and possibly two different tip voltages). 
The possibility to probe the LDoS at different spatial lo- 
cations in STM experiments offers new perspectives in 
comparison with previous experimental studies of fluc- 
tuations of the LDoS at a fixed position, using resonant 



tunneling through a localized impurity stated—. 

More explicitly, from the LDoS p(r,u), which depends 
on tip position r and voltage ui (although we keep the 
electron charge e = —\e\ and Planck's constant h in what 
follows, we assume that voltage, energy and frequency are 
loosely identified with each other), we define the nonlocal 
disorder averaged LDoS-LDoS correlations 



x(r,wi,w 2 ) = (p(ri,wi)p(r 2 ,a;2)) 



(1) 



as the centered two-point correlation function of the 
LDoS (here r = |ri — r2 1). Clearly this is a complicated 
object that depends on two tip voltages, but only on the 
distance between the two positions of the STM tip be- 
cause of translation invariance and spatial isotropy after 
averaging. 

Before turning to detailed calculations, we wish to 
give some general physical interpretation of this physical 
quantity. The basic idea is that important correlations 
are obtained whenever the two quantum cyclotron rings 
(associated to circular wave functions in a perpendicular 
magnetic field B) have a large spatial overlap (disorder 
plays, however, a crucial role in correlating locally the 
states, as we will see later on). Focusing first on the spa- 
tial dependence of Wi,^) for equal tip voltages, one 
readily understands that the intersection of the quan- 
tum rings (with width of the order of magnetic length 
Ib = \/fic/\e\B) around the two tip positions provides 

an increased area of intersection [from l B to l B (R 1 ") 1 ^ 2 ] 
when the distance between the points is close to the sum 
of the Larmor radii R L (see Fig Q] and caption for de- 
tails) . This already suggests that a peak should occur in 
the LDoS-LDoS correlations for this particular distance. 
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a) 




b) 





FIG. 1: Geometric interpretation of correlations of the LDoS 
in terms of overlap of quantum cyclotron rings of width Ib- 
Panel a) taken for close intertip distance |ri — r2 1 >C 2R L 
presents a nearly maximal overlap area of the order IbR L 
(with R L the Larmor radius), hence maximal correlations. 
Panel b) taken in the intermediate situation ri — r2 1 — R L 
has a small overlap area 1%, thus relatively weak correlations. 
Panel c) taken at tangent cyclotron orbits |ri — rg| ~ 2R L has 
a small increased overlap area of the order l 3 J 2 (R L ) 1/2 (for 
R L > Ib), so that a moderate peak as a function of intertip 
distance occurs. Finally panel d) presents the situation of 
non-overlaping cyclotron orbits at |ri — T2\ 3> 2R L , hence 
exponentially suppressed LDoS-LDoS correlations. 



model is a singlc-particlc Hamiltonian for an electron 
confined in two dimensions in the presence of both a per- 
pendicular magnetic field B and an arbitrary potential 
energy V(r), 



H 



1 



2m* 



■iV r --A(r) +V(r), 



(2) 



with the vector potential A such that VxA = B = Bz, 
and m* the electron effective mass [here r = (x, y) is the 
position of the electron in the plane] . We do not consider 
the case of graphene here, which can be easily extended 
following our previous results^, although this system is 
quite relevant experimentally for the considerations of 
the present work. 

The starting point is the realistic assumption of large 
cyclotron frequency (compared to local amplitude fluctu- 
ations of the disordered potential) at large magnetic field, 
so that Landau level mixing can be disregarded and the 
guiding center coordinate R follows a quantum motion 
along weakly curved equipotential lines. In that case the 
guiding center Green's function obeys a single pole struc- 
ture for each Landau level n as found in Refs. [ll| and [H 
[this generalizes the results of Ref. [TU to arbitrary Lan- 
dau levels]: 



We note that this effect has no classical analog, as quan- 
tum cyclotron rings collapse into circular cyclotron orbits 
at vanishing Ib (in fact, we will see that a kink instead of 
a peak occcurs in the classical limit). The energy depen- 
dence (given by the two tip voltages) can be also easily 
inferred, for instance, at coinciding tip positions. When 
both tip energies precisely match the Landau levels, max- 
imal overlap of the whole quantum cyclotron rings occur, 
leading to sharp and large positive correlations. However, 
detuning the two tip energies will probe correlations be- 
tween different circular wave functions, and destroy the 
correlations. In that case, the square averaged term in 
Eq. (P) will dominate the averaged square one, leading 
then to weak and negative correlations. 

The paper is organized as follows. Section[H]recalls the 
form of the local density of states (for fixed but smooth 
disorder) in quantum Hall systems^. We consider, in 
turn, the disorder averaged LDoS in Sec. IIII1 allowing 
us to present the technique on a simple and more usual 
quantity, and to test the approximation scheme. Then 
Sec. |TV] presents the derivation of the disorder averaged 
LDoS-LDoS correlations, which arc finally discussed in 
great detail in Sec. |V] 



II. LOCAL DENSITY OF STATES AT HIGH 
MAGNETIC FIELD 

We now present our theoretical analysis of this prob- 
lem, which can be carried out fully analytically for 
smooth disorder in the high magnetic field regime, and 
which confirms the above argumentation. The basic 



G„(R) 



UJ - E n - V n (R) + ifT 1 



(3) 



where E n — huj c (n + 1/2) are the Landau level energies 
and ui c = \e\B/(m*c) is the cyclotron frequency. One 
important quantity above is the effective potential V n (R) 
that results from averaging the bare disorder potential 
V(R) along the quantum cyclotron motion 



V n (R) = J d 2 r/F„(R 



(4) 



The kernel F n (R) here is given by the following 
expression^: 



F„(R) = 



D -A a R 2 /l 2 E 



im\l 2 B ds r 



1 



(5) 



with A s = (1 — s)/(l + s), and can also be written in an 
equivalent formic 



(-1)™ /2R 2 \ r2/ ,2 



(6) 



where L n (z) is the Laguerre polynomial of degree n. This 
kernel is also known as the form factor in the literature 
on quantum Hall effects^. Expression §5§ turns out to 
be very useful for the study of the first several Landau 
levels, while expression ([5]) is more suited for the con- 
sideration of high Landau levels (n 3> 1). Apart from 
the case n = 0, we note that the kernel F n (R,) is not a 
positive-definite function, and cannot be interpreted as a 
wave function probability density. Instead, it rather cor- 
responds to a Wigncr distribution because the physical 
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space of the guiding center coordinates R = (X, Y) is, 
in fact, associated to a pair of conjugate variables owing 
to the commutation relation [X, Y] = il 2 B in the operato- 
rial language. However, integration of the kernel F n (R) 
over an arbitrary line provides a translationally invari- 
ant Landau states probability density. In the classical 
limit n — ¥ oo, while keeping the (Larmor) cyclotron ra- 
dius = V^nlB finite (hence for Is — > 0), one getsi^ 
F„(R) ~ 2itR l *KW — Rn); 80 that the effective poten- 
tial Eq. (01 corresponds to an average over the classical 
cyclotron orbit, as previously shown in Ref. [iH . For a 
nonzero Ib, F n (R) is an oscillating function that shows a 
sharp peak of width Ib centered around |R| ~ R^. This 
quantity will be a crucial ingredient later on for the math- 
ematical identification of the quantum cyclotron rings. 

Finally, the LDoS is readily connectec&ii to the guid- 
ing center Green's function given in Eq. ([3]) 



P(r,w) 



/S|>(K-r)>^(R, 



(7) 



Here one did not include the overall spin degeneracy. In 
practice, the signal measured by local scanning tunnel- 
ing spectroscopy is directly related to the local density 
of states through an energy convolution with the deriva- 
tive of the Fermi-Dirac distribution and the experimental 
resolution window. Apart from the condition of high cy- 
clotron energy, the present calculation is valid for smooth 
potentials at the scale oils- Such a restriction on the spa- 
tial variations of the potential is, in fact, not very drastic 
because the above results are exact for arbitrary one- 
dimensional (ID) potentials^ (even very rough ones). 
For the realistic two-dimensional (2D) potential the va- 
lidity of the calculation is controlled by the small energy 
scale associated typically to potential curvature^ 



11 



urc = fV\dxxV n d YY V n - {d XY V n ) 2 \. (8) 

A more precise evaluation of the degree of reliability of 
our approximation scheme for rough potentials will be 
given at the end of the next section. Before considering 
the correlations of the LDoS, we compute now in some 
detail the disorder averaged LDoS itself. 



III. DISORDER AVERAGED LDOS AT HIGH 
MAGNETIC FIELD 

We first examine the disorder averaged LDoS, obtained 
experimentally by spatially sampling the STS current 
over a single STM tip position. Theoretically, the aver- 
aging procedure of expression ([7]) will be carried through 
an isotropic distribution function C(q) in Fourier space 
(here q = |q|) that describes the spatial correlations of 
disorder 



(y(Rx)y(R 2 )) 



C(Ri - Ra) 
d 2 q 



(9) 



(27T)' 



C(q)e 



-jq-(Ri-R 2 ) 



Typically one can take C(R) = u 2 e _ l R l 2 ^ 2 , defining 
the correlation length £ and the root mean square value 
v = \J (V(R) 2 ) of the bare disorder distribution. In that 

case C(q) = itv 2 ( 1 e-t 2q2/i . Averaging of the LDoS © 
then simply follows from exponentiating the single pole 
in Eq. ([3]) by going in the time domain, and performing 
Gaussian integration over all possible disorder realiza- 
tions 



d 2 q \V(ci)\ 2 \ 



n=0 



x exp { i[uj - E n - V n (R)]t 



(2tt) 2 C(q) 



The effective potential Eq. Q is given in Fourier space by 
V n (q) = F n (q)V (q) , where F n (q) is the Fourier transform 
of kernel ([5]) , which is easily shown to obey 



F n (q) 



d 2 re iq - r F n (r) 



1 d n 



n\ ds n (1 + s)A s 



= (-l)»7r/|f n (i|q/2). 



s=0 
(11) 



One can then readily perform the functional integral over 
the disorder realizations in Eq. ((101) : 



Mr,")) = 



x exp 



~:t 2 r 2 



(12) 



where the energy width T n is given by the relation 



d 2 a 2 
~^C(q) F n (q) 



(2tt) 2 



(13) 



This result was obtained initially in Ref. [B| for the case 
of graphene. 

Expression (fT!?]) for the averaged density of states 
(DoS) is obviously r independent, so that the R inte- 
gral can be carried using the normalization condition 
J d 2 RF„(R) = 1. The remaining time integral gives 
the final result 



(P( r ^)) 



2tt/| ^ v^r„ 

D n—0 v 



exp 




(14) 



which takes the expected Gaussian lineshape. 

The rcnormalized disorder width T n given in Eq. (|13j) 
can be analyzed^ in the classical limit n — > +oo, keeping 
the cyclotron radius R^ = ^f^nls fixed and Ib — > 0. In 
this regime, we recover results first derived in Ref.— using 
a completely unrelated method 



r„ = 



lj^2C(q)lJ a (RM 2 , 



(15) 
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FIG. 2: (color online) Disorder averaged density of states [in 
units of (2ttI%v)~ 1 ] as a function of energy for two values of 
the correlation length of the disorder distribution ({/is = 5 
and £/£b = 1) spanning the first four Landau levels. 



with Jo{z) the zeroth order Bessel function. For very 
large classical orbits such that R^ ^> £, the large- 
argument asymptotics of the zeroth order Bessel function 
Jq{z) ~ yj2/{-Kz) cos(z — 7r/4) can be used, so that 



l Jdq2C(q) 
7T 2 i?£ 



,1/4 ' 



(16) 



showing a decrease of the Landau level energy width 
with increasing index n. We note that expression (|13[) 
is more general than the classical result (fi"5f because it 
incorporates wave function spreads on the scale Ib, a 
purely quantum lengthscale which has completely disap- 
peared in the classical limit. In all cases (classical or 
quantum), the general trend is that the cyclotron mo- 
tion averages out the local potential at increasing radius 
R^, so that the energy width of the sample averaged DoS 
decreases with n. This effect is clearly seen-ir— from the 
experimentally measured spatial dispersion of the LDoS, 
which shows a rapid narrowing for higher Landau lev- 
els. However, in the opposite limit of very smooth dis- 
order £ ^ R%, this averaging by the cyclotron orbits 
becomes less efficient, and the energy width r„ depends 
very weakly on the Landau level index n. Both regimes 
are presented in Fig.f2]showing the energy-dependent dis- 
order averaged density of states for two values of £,/Ib- 

We end up by commenting on the reliability of our cal- 
culations for realistic disorders, and possibly rough ones 
at the scale of Ib- From the construction of the quan- 
tum guiding center theor y 11 ' 12 as a systematic gradient 
expansion, the lowest order calculations used here be- 
come exact for smooth disorder, namely £ Is- In 
the opposite limit £ <C Ib, the approximation scheme 
does not fully break down per se, thanks to the wave 
function averaging on the scale Ib performed within the 
effective potential (UJ) . This implies that the effective po- 
tential varies on the scale \/£, 2 + 1% and remains smooth 
even for very rough bare potential (£ <C Ib)- In this 
regime, the theory simply misses then the small param- 




-4 -2 2 4 

(to - Eq)/uj c 

FIG. 3: (color online) Comparison of the disorder averaged 
density of states [in units of (2-kI%Fo)~ 1 ] in the Landau level 
n — between Wegner's exact formula and our leading order 
gradient approximation for the extreme case of (5-correlated 
disorder (£,/Ib = 0). Rapid convergence of the approximation 
scheme is expected as soon as £ > Ib- 



cter l B /(£, 2 + 1%) °f the case £ 3> Ib, and its adequation 
becomes purely a quantitative matter. Interestingly, we 
are able to assess its validity thanks to Wegner's exact 
solution^ for the disorder averaged density of states in the 
Landau level n = 0, calculated for (5-correlated disorder 
(V(R)V(0)} =wSW(R) 



1 2 3 / 2 cxp[2(u - E o ) 2 /T 2 

2^|7r 3 / 2 r l 



{Erfi[V2(o 



£o)/r )]} 2 

(17) 



In the previous formula, Erfi(a;) = 2tt x I 2 Jjf due" is 
the so-called complex error function and Tq = w/(irl B ). 
In our calculation, the ^-correlated potential is obtained 
from the £ — > limit in Eq. ([9]) , which corresponds to the 
most stringent limit to test our approximation scheme. 
The comparison (derived for the same microscopic disor- 
der parameters) between the Gaussian expression (|14[) in 
the lowest Landau level n . = 0, using the lincwidth (|13l) . 
with Wegner's exact formula^, written in Eq. (JTTJ) , is 
given in Fig. [3] Although neither the peak value nor 
the tails (not shown) are exactly reproduced within our 
approximation, the result is close to the exact answer, 
and this is surprising because we are far from the naive 
domain of validity of the theory. 

This comparison gives some evidence that both the dis- 
order averaged LDoS and its nonlocal correlations (to be 
calculated now) , will be quantitatively obtained by a gra- 
dient expansion in the disorder for realistic cases where 
£ ^ Ib (granted by the high magnetic field limit). More- 
over, we stress that Wegner's method does not apply to 
obtain two-particle averages, and cannot be used to cal- 
culate the LDoS-LDoS correlations, showing the greater 
generality of our approach. 
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IV. DISORDER AVERAGED LDOS-LDOS 
CORRELATIONS AT HIGH MAGNETIC FIELD 

We are now ready to calculate the two-point correla- 
tions of the LDoS, taken at two different tip positions 
ri,r2 and for two different energies wi, uj 2 . The starting 
point is the expression 



'dh fdt-2 /a 2 Ri fd 2 R 2 t^ 



E 



2nl 2 J 2nl 2 

a £>n 1 =[lri2=U 

xF ni (Rx - n)F n2 (R 2 - T2 ) e K^-En^+i{^-E n2 )u 
xA(Ri,ni,ii;R 2 ,7J2,t2) (18) 
where the following disorder average must be performed: 

A(R X , m, h; R 2j n 2 , t 2 ) = fw i^+v^M 

d 2 q |f(q)| 2 l 



x cxp ■ 



(2tt) 2 C(q) 



exp 



cm 



1=1,2 



UF ni (q)e 



= exp j-I [tfr^ + t 2 2 r 2 n2 + 2tit a S niin3 (Ri - Ra)] J . 

The energy width T n was already defined in Eq. ()13[) , and 
a spatially dependent disorder correlator now appears 



(R) 



(2tt) ! 



C(g)F ni (g)F na (g)cos(q-R). (19) 



Using expressions ([5]) and (|11[) , we can perform exactly 
the integrals in Eq. (|19[) and hnd a result which can be 
written under the form 



^ni ,n 2 (-^) 



2^o 2 



ni+ri2 



exp 



(1 - si - s 2 + sis 2 ) 1 
i? 2 



e 2 + 2B Sl , S2 Z| 



,(20) 



J s 1= s 2 =0 

where we have introduced the short-hand notation 

1 - s 1 s 2 



1 - Sl - S2 



S1S2 



(21) 



Note that the disorder correlator S ni ,n 2 {R) is isotropic 
(i.e., only depends on the distance R = |R|). In addition, 
its diagonal elements taken at R = are related to the 
energy width r„ via the relation S n , n (0) — T 2 . For a 
smooth disordered potential such that £ 3> Ib, we can 
easily check from Eq. (|2H|) that the functions S ni: , l2 (R) 
are peaked at R = and decay in a Gaussian way with 
a characteristic length scale given by £ for the first few 
Landau levels. For higher Landau levels, the disorder 
correlator 5 nii „ 2 (i?) spreads over a bigger characteristic 



distance. For instance, for n% = n 2 = n 1, it has an 
extent of the order of y/ (R%) 2 + £ 2 - 

Shifting the space integrals in Eq. (fT8|) . we see that 
the quantity (p(ri, wi)/9(r 2 , u; 2 )) is clearly described by 
a function of r = ri — r 2 only. In the following we will 
consider the centered two-point correlator of the LDoS 
defined in Eq. (JlJ. Integration over the time variables 
in Eq. (|18p can be performed analytically and yields the 
expression for the LD0S-LD0S correlator 



X(r,wi,w 2 ) 



2tt1 2 J2tt1 2 ^ ^ 



F ni (Ri-n) 



xF„ 2 (R 2 -r 2 )9 ni ,„ 2 (|Ri-R 2 |, Wl ,w 2 ) (22) 



with 



e nii „ 2 (i?,wi,o; 2 ) = - [[r^r 



,(R)H 



[■Tni 1^712] e 



}. (23) 



Although the above expressions are still too complicated 
to make precise statements on the nature of the LDoS- 
LDoS correlations, we can already infer some of the early 
predictions formulated in the Introduction. Clearly the 
two kernels F n (R — r) in Eq. (|2"2"|) centered on the two po- 
sitions of the tip ri and r 2 impose a constraint on the two 
guiding center coordinates Ri and R 2 , which must live 
predominantly within cyclotron rings of extent R^ and 
i?^ 2 with width Ib (we note that the kernels have rather 
n oscillations within a given radius R^, so that there 
are in total n different rings for a given quantum state). 
This formula thus confirms our initial expectation that 
the area of overlap between two cyclotron rings dictates 
the behavior of the LD0S-LD0S correlations. A second 
interesting aspect is that it is the additional disorder- 
dependent kernel Q ni ,n 2 which permits nontrivial corre- 
lations. Should this complicated function of energy and 
space be negligible, one would end up with vanishing cor- 
relations. 

Yet, to get more quantitative understanding of the 
complete integral in Eq. (|22[) . one needs to further sim- 
plify the above expressions. This can be achieved by in- 
troducing the change in variables R' = (Ri + R 2 )/2 and 
R = R 2 — Ri . After integration over the center-of-mass 
position R', we get 

X(r,wi,w 2 )= E /^-7 2 -/™i 1 «2( R + 1 ') 

n 1= n 2 =0 J 77 B 

xe ni , n3 (fl,m,w 2 ) (24) 



with 



/ni,n 2 (R-) -/ 27rf 5 ! F «i 



R 



R 



F nn R' 



R 

~2~ 



(2Q- 2 
ni \n 2 l Os^ds™ 2 



d n 1+ n 2 C Xp(-R 2 /[2^S Sl , S2 ]) 



1 - SlS 2 



(25) 



si=s 2 =0 
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where we have used expression ((5j) for the kernel F„(R). 
Then, introducing the polar coordinates for the position 
R and performing the angular integral in Eq. (|24|) , we ar- 
rive at the final expression of the disorder averaged LDoS- 
LDoS correlations (which now obviously is only function 
of the intertip distance r = |r|): 

-f OG +00 



Eq. (|23|) . Therefore, the main contributions to the cor- 



r+00 

X(r,wi,w 2 )= "ll l B 2 RdRh nun2 (R,r) 

n 1= 0n 2 =0 ^° 

(i?,wi,<j 2 ), (26) 

where 

h nun2 (R,r) = J — /„ 1>n2 (V-R 2 + r 2 + 2rRcos 9 J 
_ (2tt/|)- 2 d ni+n ' 2 



ni ln 2 l ds^ds™ 2 

exp(- [R 2 +r 2 ]/[2P B B SuS2 ]) 
1 - sis 2 



rR 

B Si >S2 



(27) 



si=s 2 =0 



with /o(^) the modified Bcssel function of the first kind. 
An alternative writing of Eq. ([2"7| obtained by going to 
the Fourier space, which turns out to be more suitable 
for the consideration of large «i and n 2 , is as follows: 



;2 

l B 



(2wi%y Jo 



dq qF ni (q)F n2 (q)J (Rq) 

x J Q {rq). (28) 



The above expressions (|26|) through (|28p are exact for 
arbitrary 2D smooth disorder (i.e., £ 3> is). 



V. DISCUSSION OF THE LDOS-LDOS 
CORRELATIONS 

A. Spatial dependence of the correlations 

Let us now analyze the LD0S-LD0S correlations on 
the basis of Eq. (j2l))) . We see that the correlations result 
from the combination of two functions h ni ^„ 2 (R,r) and 
,712 (R,u)i,u>2), which are quite different in nature. 
Obviously, only the functions 6 nii „ 2 (R,uji,ui2) defined 
in Eq. (|23p contain the information about the energy de- 
pendence. For sufficiently widely separated Landau lev- 
els, the overlap between two states with two given ener- 
gies u>i and L02 vanishes for most of the Landau level pairs 
due to the sharp energy Gaussians cutoff in Eq. (j2U|) , so 
that essentially only the one specific pair of Landau lev- 
els (711,712) that is associated to the cyclotron energies 
(E ni ,E n2 ) closest to the tip voltages (wi,^) yields a 
nonzero contribution in the sum over Landau levels in- 
dices in Eq. ([26]) . Furthermore, the LD0S-LD0S corre- 
lator x(r, wi,W2) clearly vanishes whenever the disorder 
correlator S nitn2 (R) is small compared to r ni r n2 due to 
the exact cancellation by the square averaged term in 



relations arise when S„ 



>.(R) 



r r 

x ni x TI2 ■ 



In contrast, the functions h ni . n2 (R,r) are indepen- 
dent of the characteristic features of the disorder and 
have a pure geometric origin (as already discussed, they 
contain information on the spatial overlap of the quan- 
tum cyclotron rings). For r = 0, we have the relation 
h ni ,n 2 (Ri 0) = / nii „ 2 (R) from which a simple physical 
interpretation can be easily drawn. Using the semiclas- 
sical physical picture of the kernel F n (Il) put forward 
previously, we understand that nonzero contributions for 
/ nii „ 2 (R) in Eq. (|25[) arc picked up only from the (quan- 
tum broadened on scale Ib) cyclotron orbits with radii 
R^ and R^ 2 that live around the points ±R/2 sepa- 
rated by the distance R. This statement can be proved 
on more mathematical grounds starting from expression 
([251) . Taking the limits m and n 2 — > 00 in Eq. ([2"g)l 
as done in Ref. [TBI , we obtain an approximate formula, 
which reads 



h nu n 2 (R, r) 



72 

1 B 



(2^r ^0 



dqqJoiR^q) J (R^ 2 q) 

xJ Q (Rq)J (rq). (29) 



For r = 0, integral ([29)) strictly vanishes when R > 

Rn, + Rt_ or R < 
istic case of finite 



\Rn! 



R^., I ■ In the more real- 



fii and ri2, this semiclassical limit 
shows that the resulting overlap of two quantum orbital 
motions turns out to be significant under the inequal- 



\R„ 



Ri 



< R < R r 



Ri 



The functions 



7ii.7i2 (R-i r ) being symmetrical in R and r, we understand 



lties 
h 

that for R = the contributions to the LD0S-LD0S cor- 
relations arise when the distance r between the two tip 
positions is such that | R% 1 — R^ | < r < R^ i + R^ 2 , 
precisely as anticipated in the discussion of Fig. [TJ 

The spatial dependence of the two-point correlator 
x(r,LUi,U2) resulting from the numerical computation of 
the integral over the distance R in Eq. ([26]) is shown 
in Fig. 2] in the regime of widely separated Landau 
peaks (we have taken in what follows hw c = 5v). In 
these two figures, we have chosen the situation where the 
LD0S-LD0S correlations are maximal, that is, we have 
considered that the energies u)i (here i = 1,2) corre- 
spond exactly to Landau level energies E nt . The case 
of equal energies u>\ = lu 2 = lo is first investigated in Fig. 
21 According to the previous discussion, the dominant 
contributions among the different possible pairs of Lan- 
dau level indices are the diagonal ones corresponding to 
m = jt.2 = n. For the lowest Landau level energy (n = 0) 
shown with the solid line in Fig. 21 the correlations de- 
crease in a monotonic way as a function of the distance r 
between the two tip positions with a characteristic decay 
length of the order of the disorder correlation length £ 
(here all the distances are expressed in units of the mag- 
netic length Ib and we have taken £ = 51b)- When the 
energy of the first Landau level is probed (dotted line of 
Fig. |H corresponding to u = E\), the spatial dependence 
of the correlations is still decreasing with the same decay 
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length £, but it now exhibits a mild peak for the position 
r close to 2R^ = 2^/2uIb ~ 2.81b for n = 1, associated to 
the overlap of the quantum cyclotron rings for the n = 1 
Landau states. In the second Landau level (at uj = E2), 
besides the peak close to 2R^ = 41b for n = 2, an ad- 
ditional peak is clearly seen in the r dependence of the 
LD0S-LD0S correlations, see dashed line in Fig. SJ This 
is because the wave functions for n > 1 have n zeros, 
hence additional rings of high probability density. 



1 ' 1 ' r 




FIG. 4: (color online) LD0S-LD0S correlations at equal en- 
ergies (cji = U12 = w) as a function of the tips distance r 
(in units of Ib) for different energies u>. Here the correla- 
tions are expressed in units of (2irl%v)~ 2 with the parameters 
Tvuj c /v — 5 and Zb/£ = 0.2. The spatial dependence of the 
correlations presents peaks whose number corresponds to the 
Landau level index of the probed energy (for instance, for 
uj — E2 two peaks are seen). 



This nonmonotonic spatial dependence observed for 
equal energies was anticipated in the Introduction of the 
paper, and can be understood more quantitatively as fol- 
lows. When the disorder potential is smooth on the scale 
of Ib (i.e., £ 3> Ib) the two functions involved in the 
integrand of Eq. ((25)) are, in fact, characterized by two 
very different characteristic length scales. Indeed, the 
function h ni ^ n2 (R, r) decays with R on the typical length 
scale R^ = y/2nl b which, for indices n not too big, turns 
out to be much smaller than the characteristic length 
scale (of the order of \J (R^) 2 + £ 2 ) for the spatial vari- 
ations of the functions ® n ,n {R, w, w). Therefore, within 
the first few Landau levels a good approximation to the 
integral ([26]) is provided by the Laplace's method owing 
to the inequality £ 3> R^. Using the formula 




we obtain the approximate analytical expression for the 
LD0S-LD0S correlations at identical tip energies (oj\ = 
L02 = uj) for uj close to the energy E n 



x(t,uj,uj) 



(2irl 



2irv 2 



l V^C e -( u -E n f/(2v*) 
\2V21b 



& 1 



(n!) 2 ds^dt 



1 - sis 2 



zL/2 R 

B s l j s 2 



e 

Si=S 2 =0 



(31) 



Consequently the nonmonotonous behavior of the cor- 
relations seen in Fig. |4] is connected with the oscilla- 
tions of the function h n>n (0,r) within Eq. (|2l)]) . and is 
fully described analytically by formula (PTj) in the regime 
£ 3> Ib- As seen in Fig. [5] for two large values of £/Zb and 
for energies taken at the second Landau level, the agree- 
ment between the numerical evaluation of the integral in 
Eq. (|26p and the analytical approximation pip is quite 
excellent. We note also that the correlations disappear 
(for generic values of the frequency) in the clean limit 
v — > despite the v~ 2 prefactor due to the exponential 
terms in Eq. (|3I[) . while they become singular precisely 
at the Landau level frequency. 

We would like to comment here on the semiclassical 
limit of fixed R^ with n — > 00 (and vanishing Ib) for the 
LD0S-LD0S correlations. In that case, we replace the 
kernel F n (R) by ^-L^flRl-i^) in Eq. (JUJ), and obtain 



numerically Fig. [6j Not surprisingly, the quantum oscil- 
lations at r < 2R^ disappear, yet a kink subsists for the 
exactly tangent cyclotron trajectories at r = 2R^, while 
a logarithmically diverging correlation occurs at r = 0. 
These are classical remnants of the quantum effects dis- 
cussed so far, and this robustness can again be expected 
according to the geometrical argument of Fig. [T] (now 
in the classical limit). These characteristic features can 
be understood analytically. Indeed, inserting Eq. (|29|) 
within formula (|26p . and performing Laplace's method 
on the independent integral (this is valid in the case 
£ > R L ), we find 

X (r,E n ,E n ) ~ ^^l(iJ%J (rq)[J (R L q)] 2 -l\ 

(32) 

Expression (|32p is clearly logarithmically divergent at 
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FIG. 5: (color online) LDoS-LDoS correlations at equal en- 
ergies taken at the second Landau level (lui — 0J2 — E2) as 
a function of the tips distance r (in units of Ib) for two dif- 
ferent values of the disorder correlation length £/£b = 14 and 
£/Ib = 5, with a comparison of the numerical evaluation of 
expression (|26[) and the analytical formula (|31[) . 




r/K 



FIG. 6: (color online) LDoS-LDoS correlations at equal ener- 
gies taken at the Landau level (ui\ = 102 = E n ) as a function 
of the tips distance r (in units of R^) in the semiclassical 
approximation (Ib — > 0) to integral (|22|l . 



Ea we get X (0,E n ,E n ) oc log(n)/i?£ with R% = y/2nl B . 
We recover the logarithmic divergence for n — > 00 and 
fixed R~n, while this expression vanishes as \og(n)/^/n 
for 71 — > 00 and fixed , in agreement with Fig. Ul where 
a slight decrease with increasing n is seen at r = 0. We 
thus stress that quantum mechanics (finite Is) always 
regularizes the singular classical behavior, and that max- 
imal correlations are, in fact, obtained (at a given mag- 
netic field, hence fixed Ib) for the lowest Landau levels. 

Let us come back to the quantum case of finite I b and 
investigate the effect of energy detuning in the spatial de- 
pendence of the LDoS correlator. For energies lo\ 7^ 
a different spatial structure from the situation of equal 
energies can be seen in Fig. which corresponds to 
the cases (wi,w 2 ) = (E ,Ei) and (cJi,w 2 ) = (Ei,E 2 ), 
the solid and dotted lines, respectively. Unlike the equal 
energy cases of Fig. [5J the correlations arc no more max- 
imally obtained for the tips distance r = but for an 
intermediate tip distance of the order \R± — R%\, as can 
be guessed again from the geometrical interpretation of 
Fig. [1] in the case i?f 7^ R 2 ■ We note that the spatial 
dependence of the correlations between the first and sec- 
ond Landau levels is also characterized by an extra mild 
peak, a reminiscent feature of the equal energy case. 




r = 0, but continuous at r = 2R L (the situation of 
tangent cyclotron orbits). The spatial derivative of the 
correlator can be analyzed in the limit r —> 2R L , and 
gives a finite derivative for r < 2R L and a diverging one 
for r > 2R L , explaining the kink feature seen in Fig. [6] 
The onset of the logarithmic divergence at r = can be 
understood starting from the quantum expression pip , 
which alternatively reads for r = and u> = E n : 



X (0,E n ,E n ) 



(2nll) 



2irv 2 



V21 



BJO 



du \L r , 



(33) 



The above equation is again obtained from Laplace's 
method in the case £ 3> -R„ to analytically perform the 
integral over R in Eq. (|22|) , and used the explicit expres- 
sion ([6]) for the kernel F n , instead of the derivative trick 
performed in Eq. pip . Using the asymptotic formula for 
the Laguerre polynomial in the large n limit as in Ref. 



FIG. 7: (color online) LDoS-LDoS correlations at unequal 
energies (wi = uj and u?2 — u + u) c ) as a function of the tips 
distance r (in units of Ib) for different energies ui. We used 
the same parameters as in Fig. [4] 



B. Energy dependence of the correlations 

We now study the energy dependence of the LDoS- 
LDoS correlations. We have represented in Fig. [8] the 
correlator x(r, wi, W2) as a function of the energy L02 = w 
for identical tip positions (i.e., for r = 0) when the first 
tip energy is pinned to the first Landau level (lo\ = E\). 
We note the presence of different peaks in the correlations 
corresponding to the sequence u> = E n of the Landau lev- 
els energy peaks in the local density of states. However, 
it is worth noting that the LDoS-LDoS correlation peaks 
are much sharper than the LDoS peaks because of the 
sensitive matching of the spatial overlaps between the 



9 




To obtain this expression, we have developed the dis- 
order correlator S ni . ri2 (R) [Eq. (fT9|) ] and the linewidth 



2 

uj/uj c 



FIG. 8: (color online) LDoS-LDoS correlations for the tip dis- 
tance r = and for the energy wi pinned to the first Landau 
level (oil = Ei) as a function of the energy 102 = The other 
parameters are the same as in Fig. [4] The highest peak is 
obtained for loi = uji and is accompanied by satellite peaks 
corresponding to the neighbouring energy levels of the first 
Landau level. 



r„ in the large £ limit. As can be seen in Eq. (|34|) . a 
competition between two exponential terms occurs de- 
pending on the relative values of ui\ and ui2- For Suj = 
lji — LU2 + E n2 — E ni = 0, the first term in the right-hand 
side (r.h.s) of Eq. ([34|). which is large and positive (of 
the order of £/-R), dominates over the second exponential 
contribution, and one gets positive LDoS-LDoS correla- 
tions, as already noted. For 5ui slightly different from 
zero, the first term decreases exponentially fast (note the 
large £ 2 /R 2 prefactor within the exponential), resulting 
in the sharp peak observed in Fig. [8] near coinciding en- 
ergies. At increasing \5u>\, the second term in the r.h.s of 
Eq. (|34[) . which has a small and negative amplitude, is 
characterized by a slower exponential decrease, and then 
gives the main contribution to the correlations. Coming 
back to the initial Eq. ([1} for the correlator x(r, u>i, ui^), 
we conclude that the square average density dominates 
quickly at increasing energy detuning. 



quantum rings associated to the states in Landau level 
n\ = 1 and n 2 = n - Besides the strongest peak ob- 
tained for ui = Ei, a succession of lateral peaks occurs 
for all other cyclotron energies E n , showing that nondiag- 
onal elements of h nin2 (R, r) are also strongly correlated. 
More interestingly, the regions surrounding each of these 
sharp peaks with positive correlations are characterized 
by negative values for the correlator x(0, Ei,ui2), thus 
corresponding to anticorrclations. This crossover from 
positive to negative correlations as a function of the en- 
ergy can be easily understood on the basis of Eq. (|23|) at 
£ ^S> R L , which is then well approximated by a small- R 
expansion 



(i?,wi,o;2) 



1 / £ £ ("l-^2 + -E n2 -E n i> 2 



-e 4h 2 



xe 



3^ 



2irv 2 \ 2R 

(^ 1 -E„ 1 ) 2 +(u 2 -E„ 2 ) 2 



(34) 



Analytical insight can be obtained by further assum- 
ing the classical regime of R L Ib (hence high Landau 
levels). Using then the approximation (|29|) within for- 
mula (f2l))) together with Eq. (fM)) (which is valid in the 
limit £ 3> R L ) we get the semiclassical approximation for 
the LDoS correlations at coinciding tip position r = 



X(0,wi, uj 2 ) 



2nv 2 {2irl 



-(p(0,uji))(p(0,u} 2 )). 



J2 / dqqJ (R^q)J (R^q) 



dRJ (Rq)e~ 



l)("2- 



(35) 



Here the double sum over (711,712) is again constrained 
by the external frequencies (u)i,u>2), and reduces to a 
single term in case of sharply defined LDoS-LDoS corre- 
lation peaks, giving rise to a single peak at uji — w 2 = 
E ni — E n2 . Let us focus first on the case where ni = n 2 , 
so that = Rfl 2 = R L . We see that the above in- 
tegrand in Eq. (|35[) then behaves as 1/q for vanishing 
Sui, because \Jo{R L q)\ 2 looses its oscillatory character 



at large momentum q. One obtains thus a logarithmi- 
cally diverging peak x(0, 1^1,1^2) ~ log |«/5a;|, related 
to the log \Ib/t\ spatial divergence found previously in 
the semiclassical limit at small intertip distance. How- 
ever, for nonzero n = ri\ — n 2 , the LDoS-LDoS correla- 
tion peak at u>i — uj 2 = E n is no more logarithmically 
diverging when Sui — > 0. This is because the product 
Jo{Rn 1 q)Jo(Rn 2 q) behaves as cos(l 2 3 nq/R L )/q for large 
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momentum q, as can be seen using the expression for the 
cyclotron radii R^ ± 2 = ^^^r in the limit n <C n\, n 2 , 
where R L = \jn\ + ^Ib- In order to capture the rele- 
vant energy scale at small Suj, we can make the change 
in variables k = (^\Suj\/v)q and z — (v/£,fl n )R with 
f2„ = nvl 2 B /(R L £) into Eq. and use the asymptotic 
form of the Bessel functions for \5u\ <C £v/R L . After in- 
tegration over k, this provides the deviation of the LDoS 
correlations from the n th peak value (at Slj = 0) 



6 X 



(2^|)- 2 g 
2tt 2 v 2 R l 



+00 



dz 



1 



(36) 



The new energy scale f2„ therefore sets the width of 
the n th correlation peak. Since f2„ <C v in the regime 
we recover the fact that the LD0S-LD0S 
correlations are more sharply defined than the average 
LDoS peaks. The linear increase of tt n with n also ex- 
plains the progressive smearing of the correlations at in- 
creasing energy detuning of the tips (sec Fig. [5J. 

We now discuss the situation of strongly overlapping 
Landau level, cu c <C v = y/ (V 2 ) , so that the average den- 
sity of states becomes structureless. In that case, the 
sharper peaks in the LD0S-LD0S correlations can sur- 
vive for smooth disorder (£ ^> lg), under the condition 
vIb/£, *C w c for the lowest Landau levels. A similar re- 
sult was already found by Rudin et al^, who studied the 
energy dependence of the local LD0S-LD0S correlations 
in a weak magnetic field and long classical orbits R L ^> £ 
in the diffusive regime. The disorder dependence of cor- 
relator (|3"Tj) in Ref. [I?] was also in 1/v 2 as in the present 
paper. 

It is worth stressing that the condition to obtain sharp 
peaks in the LD0S-LD0S correlations corresponds pre- 
cisely to the absence of local Landau level mixing. In- 
deed, transitions between adjacent Landau levels pro- 
vide^ a typical energy scale l B \ W\ 2 /oj c , so that our cal- 
culation is controlled when the parameter (l B /^ 2 )(v 2 /lu 2 ) 
is small. Clearly, large overlap in the average DoS can 
be compatible with no mixing, because this quantity re- 
lates to global properties (long-wavelength fluctuations) 
of the smooth disorder. In contrast, only the correlations 
of the LDoS can feel the local interplay of disorder and 
Landau quantization, and reveal whether the Landau in- 
dex n stays a good quantum number or not. 

Interestingly, a logarithmic singularity for the n = 
peak and an energy width f2„ proportional to n for 
the n th peak were obtained in the semiclassical diffusive 
regime- 7 . This is only slightly different from our results 
in the semiclassical limit, showing the continuity of the 
present physics from low to high magnetic fields. We 
emphasize again that quantum effects at finite magnetic 
length regularize the spurious divergences generically re- 
lated to the semiclassical approximation. 

Finally, we note that some of the energy-related fea- 
tures on the LDoS correlations discussed above have al- 
ready been reported experimentally^ - in heavily doped 
three-dimensional (3D) GaAs semiconductors. For in- 



stance, conductance anticorrelations were observed^ - in 
the presence of Landau levels, while the narrowing of 
the LDoS fluctuations has been found^ at weak mag- 
netic field. However, these experimental studies were 
performed using resonant tunneling impurity, hence at 
a fixed position, asking for a different analysis from what 
was performed here (besides the 3D character of the stud- 
ied samples). Indeed, averaging was performed either 
over the applied voltage^ or magnetic fieldiS, and in the 
later case, the Landau levels are simply lost. Spatial av- 
eraging in a STM configuration, as proposed in our work, 
should allow a greater control of the LDoS correlations 
(in terms of the applied magnetic field and tip voltages) . 
Moreover, this offers a way to investigate spatial cor- 
relations of the LDoS, that could not be assessed with 
previous experimental techniques. 



VI. CONCLUSION 

We have studied theoretically the two-point correla- 
tions of the local density of states in a disordered two- 
dimensional electron gas under a large magnetic field. 
A rich spatial dependence of the correlations was found, 
which can be qualitatively explained by geometrical over- 
laps of the two quantum cyclotron rings that roughly de- 
scribe circular wave functions. The energy behavior of 
the correlations was shown to provide sharp peaks when 
the frequency detuning matches integer multiple of the 
cyclotron frequency, similar to the low magnetic field re- 
sults of Rudin et aliH. These sharp peaks in the LDoS 
correlations reveal that Landau levels correspond to well- 
defined quantum numbers, information that cannot be 
easily gathered from the average density of states only, 
where large overlaps are usually reported experimentally^ 
even at large magnetic fields. We have also emphasized 
here that LDoS correlations can be either positive or neg- 
ative, depending on the degree of frequency mismatch. 
We would also like to mention that at energies strictly co- 
inciding with the centers of Landau levels (quantum Hall 
critical point), the LDoS correlations exhibit a power-law 
behavior at long distance, reflecting the multifractality of 
critical wave function s 19 ' 20 . This more complex behav- 
ior is beyond the scope of the present paper, which did 
not consider quantum tunneling between closed drift or- 
bits. Such effects are, however, only relevant at very low 
temperatures for a smooth potential. 

We end up by noting that recent experimental pro- 
gresses for electron gases confined at the surface of InSb 
semiconductor— and also in graphene^ allow high spatial 
and energy resolution measurements of the local density 
of states. The disorder averaged LDoS and its corre- 
lations can, in principle, be straightforwardly obtained 
from the experimental data by averaging over large scale 
spatial maps. Because both the width of the disorder dis- 
tribution and the correlation length of the random land- 
scape can be extracted from the knowledge of the average 
LDoS, our predictions could be tested even quantitatively 
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without extra fitting parameter. In the case of graphene, 
extension of the present work can be straightforwardly 
done following the results of Ref. H). The spinorial form 
of the wavefunction is expected to lead to more complex 
spatial structures, yet with similar general behavior to 
that discussed here. 
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